Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

The entropy of Garfinkle-Horne dilaton black hole due to arbitrary spin fields

Using the membrane model which is based on brick wall model, we calculated the free energy and entropy of Garfinkle-Horne dilatonic black hole due to arbitrary spin fields. The result shows that the entropy of scalar field and the entropy of Fermionic field have similar formulas. There is only a coefficient between them.     

高折射率椭圆芯布拉格光纤的偏振特性

应用带完全匹配层边界条件的全矢量有限元方法,分析了高折射率椭圆芯布拉格光纤的偏振特性,详细讨论了光纤结构参量对模式双折射度以及群速度走离的影响,结果发现:高折射率椭圆芯布拉格光纤的模式双折射度可达10-2量级,比传统保偏光纤至少高出一个量级,并且表现出不同于传统保偏光纤的群速度走离特性;在保持较高双折射度的同时,通过合理的结构设计,可在一定的波长处灵活地获得较大的群速度走离或零走离特性,具有重要的潜在应用价值;最后,简要分析了低折射率区域的折射率变化对偏振特性的影响;研究结果有助于设计高性能的微结构保偏光纤。     

CFD在燃煤细粒子凝聚过程中的应用

燃煤产生的细粒子富集了大量的有毒痕量元素,对大气环境和人类健康造成严重危害。本文针对燃煤细粒子的形成过程,在计算流体力学(CFD)的基础上,结合气溶胶动力学理论,模拟了烟气中细粒子在圆柱形流场中由于碰撞而凝聚的过程,并通过计算结果分析细粒子颗粒特征参量(速度、质量、直径和颗粒数目)随流场的变化趋势,为深入研究燃煤细粒子的形成演化机理提供理论基础。     

Numerical Prediction of the Strain State of Thermoset Matrix/Mineral Filler Composite Plates

The numerical prediction of the fields of inelastic strains (the linear invariant of the tensor of inelastic strains) in thermoset polyester/marble filler composite plates is discussed. A uniformly distributed load is applied to the plates, which lie on a steel base. The strain fields are predicted by means of the boundary element method by using creep test data for the composites and the polyester matrix itself. Identical creep tests were performed for two ages of the materials (1 month and 13 years), which allowed evaluating the aging effect. The study is carried out in two stages. At the first stage, the application of the generalized Maxwell-Gurevich equation to the thermoset matrix/mineral filler composite is demonstrated. The model parameters determined from the experimental creep data is used for the second stage, where the state of inelastic strains in the plates is predicted by applying the boundary element method. The influence of composite formulation (filler content) and physical aging of the polyester matrix on the state of inelastic strains in the plates is shown.__________Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 41, No. 2, pp. 145–156, March–April, 2005.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

能动反射镜低频像差补偿能力分析

 采用Algor软件包并利用有限元法对37单元能动反射镜低频误差补偿能力进行了分析计算，确定了镜面厚度8mm，驱动器间距46mm，极头直径12mm等能动反射镜制造的关键参数。通过计算机仿真分析证明，基于这些参数所设计的能动反射镜可以对低频误差进行有效校正，校正精度的均方根值小于等于λ／８，从而可以为能动反射镜的制造提供理论依据。     

Analysis of end point vibrations of a two-link manipulator by integrated CAD/CAE procedures

In this work, the effect of flexibility on the trajectory of a planar two-link manipulator is studied using integrated computer-aided design/analysis (CAD/CAE) procedures. The solid models and finite element models of the parts of the manipulator are created by using the CAD/CAE software I-DEAS. The assembly is defined, and knowing the payload and the end point trajectory, the velocities and accelerations of the parts, joint forces and driving torques are calculated using the rigid body dynamics. All the time dependent nodal forces acting on the parts including distributed gravity and inertia forces are created in files with the I-DEAS program file format. The finite element vibration analysis of the parts is performed by I-DEAS. The end point vibrations and the deviations from the rigid-body trajectory are analyzed for different types of end point acceleration curves. A circular trajectory is considered as an example. It is observed that the precision of the manipulator can be increased by testing different end point acceleration curves without changing the trajectory and the duration of the end point work. The procedure explained in this work can be used for this purpose successfully.     

Sufficient conditions of non‐uniqueness for the Coulomb friction problem

We consider the Signorini problem with Coulomb friction in elasticity. Sufficient conditions of non‐uniqueness are obtained for the continuous model. These conditions are linked to the existence of real eigenvalues of an operator in a Hilbert space. We prove that, under appropriate conditions, real eigenvalues exist for a non‐local Coulomb friction model. Finite element approximation of the eigenvalue problem is considered and numerical experiments are performed. Copyright © 2003 John Wiley & Sons, Ltd.     

美式期权定价中非局部问题的有限元方法

在本文中 ,我们关心的是美式期权的有限元方法 .首先 ,根据 [4 ]我们对所讨论的问题引进一个新奇的实用的方法 ,它涉及到对原问题重新形成准确的数学公式 ,使得数值解的计算可以在非常小的区域上进行 ,从而该算法计算速度快精度高 .进而 ,我们利用超逼近分析技术得到了有限元解关于 L2 -模的最优估计 .